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ABSTRACT 

We present a theoretical investigation of CN and HCN molecule formation 
in dense interstellar clouds. We study the gas-phase CN and HCN production 
efficiencies from the outer photon-dominated regions (PDRs) into the opaque 
cosmic-ray dominated cores. We calculate the equilibrium densities of CN and 
HCN, and of the associated species C"*", C, and CO, as functions of the far- 
ultraviolet (FUV) optical depth. We consider isothermal gas at 50 K, with 
hydrogen particle densities from 10^ to 10^ cm~'^. We study clouds that are 
exposed to FUV fields with intensities (at 1000 A) from 5 x 10~^^ to 5 x 10~^^ 
erg s~^ cm~^ Hz sr~^, or 20 to 2 X 10^ times the mean interstellar FUV inten- 
sity. We assume cosmic-ray H2 ionization rates ranging from 5 x 10^^^ s^^, to an 
enhanced value of 5 x 10^^^ s^^. We also examine the sensitivity of the density 
profiles to the gas-phase sulfur abundance. 

Subject headings: galaxies:ISM - ISM:evolution - molecular processes 

1. Introduction 

Millimeter-wave line emissions of CN and HCN molecules are widely used probes of dense 
molecular gas and of photon-dominated regions (PDRs) in the Galactic interstellar medium 
(ISM). Observed sources include molecular interfaces in star-forming regions (Greaves & 
Church 1996; Simon et al. 1997; Young Owl et al. 2000; Savage et al. 2002; Schneider et 
al. 2003; Johnstone et al. 2003), refiection nebulae (Fuente et al. 1993, 1995, 2003; Jansen et 
al. 1995), planetary nebulae (Bachiller et al. 1997a), and circumstellar envelopes and disks 
(Wooten et al. 1982; Truong-Bach et al. 1987; Bachiller et al. 1997b; Lindqvist et al. 2000; 
van Zadelhoff et al. 2003; Thi et al. 2004). Observations of CN and HCN have recently been 
used to study PDRs in the starburst galaxy M82 and in other external systems (Fuente et 
al. 2005; Meier & Turner 2005). 
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In nearby objects, such as the Orion Bar and the reflection nebulae NGC 2023 and 7023, 
the molecular emissions have been mapped across the PDRs. The CN/HCN intensity ratios 
are largest near the stellar sources of the illuminating far-ultraviolet (FUV) radiation fields, 
and the intensity and density ratios decrease with increasing optical depth and distance from 
the stars. This behavior is broadly consistent with theoretical expectations (Sternberg & 
Dalgarno 1995 [SD95]; Jansen et al. 1995) and is evidence of selective photodissociation of 
HON versus CN (van Zadelhoff et al. 2003; Thi et al. 2004). 

Many discussions of gas-phase nitrogen chemistry in molecular clouds have been pre- 
sented in the literature (e.g., Hcrbst & Klcmperer 1973; Prasad & Huntress 1980; Viala 
1986; Pincau dcs Forets et al. 1990; Herbst et al. 1994; SD95; Lee et al. 1996; Turner et 
al. 1997). Comprehensive reviews of PDR observations and theory, and related subjects have 
been presented by HoUenbach & Tielens (1999), Sternberg (2004), and van Dishoeck (2005). 

In this paper we focus on the gas-phase production of CN and HCN, and present results 
for a wide range of conditions. We analyze how the CN and HCN formation and destruction 
sequences vary with optical depth, flrst through the PDRs and then into the opaque cosmic- 
ray dominated cores. We discuss how the (equilibrium) density profiles depend on the cloud 
hydrogen gas densities and the incident FUV field intensities, and we identify the qualitative 
changes in the density profiles that may be expected in moving from low- to high-density 
systems. The production efficiencies of the carbon bearing CN and HCN molecules depend 
on the availability of free C^ ions and C atoms in the gas, and we present and discuss the 
associated C"*", C, and CO density profiles in the parameter space we consider. In this paper 
we do not consider HNC or the isomeric abundance ratio HCN/HNC (e.g., Watson 1974; 
Schilke et. al. 1992; Herbst et al. 2000) which, together with CN/HCN, may be expected to 
vary with optical depth, density, and FUV field strength. 

In §2 we describe the basic ingredients of our models. In §3 we discuss the gas-phase 
CN and HCN reaction sequences that operate in molecular clouds. In §4 we present detailed 
results for a "reference model" to illustrate the depth-dependent formation pathways and 
density profiles. In §5 we present our parameter study and discuss results for a range of 
gas densities, FUV field strengths, cosmic-ray ionization rates, and gas-phase elemental 
abundances. We summarize in §6. 

2. Model Ingredients 

We performed our model computations using an updated version of the SD95 code. The 
models consist of static, plane-parallel, semi-infinite slabs, exposed on one side to isotropic 
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FUV (6-13.6 eV) radiation fields. The steady-state abundances of the atomic and molecular 
species are computed as functions of the visual extinction, Ay, from the cloud surface. 
The models account for scattering and absorption of the FUV photons by dust grains, and 
the resulting depth-dependent attenuation of the atomic and molecular photodissociation 
and photoionization rates. The effects of H2 and CO absorption-line shielding arc also 
included. The molecular chemistry is driven by the combined action of FUV photoionization 
and photodissociation, and cosmic-ray impact ionization. The resulting sequences of (two- 
body) gas-phase ion-molecule and neutral-neutral reactions are mediated by dissociative 
recombination and photo-destruction. 

We assume that the spectral shapes of the incident FUV fields arc identical to the Draine 
(1978) representation of the interstellar field in the solar neighborhood (see also Parravano, 
HoUenbach & McKee 2003). The field intensity is parameterized by a scahng factor x, where 
for X = 1 the FUV intensity at 1000 A is 5.4 x 10~^° erg s~^ cm~^ Hz~^ sr~^. We construct 
models for x ranging from 20 to 2 x 10^, appropriate for PDRs in the vicinity of young OB 
stars and clusters (Sternberg et al. 2003). 

In our computations we include the same set of 70 atomic and molecular carbon, nitro- 
gen, oxygen, sulfur, and silicon bearing species ^ considered by SD95. Complex hydrocarbons 
and PAHs are excluded. We solve the depth dependent equations of chemical equihbrium 

Yl hji{T)njni + ^[x^ij + iij\nj = |^ kjuui + ^[x^ji + iji] } (1) 
jl j jl j 

as functions of Ay. In these equations, are the densities (cm~^) of species i, and kiji{T) 
are the (temperature-dependent) rate coefficients (cm"^ s^^) for chemical reactions between 
species j and / that lead to the production of i. The parameters x^ij ^^'^ iij ^^^^ the 
FUV photon and cosmic-ray destruction rates (s~^) of species j with products i. We solve 
Equations (1) via Newton-Raphson iteration, and use Bulirsch-Stoer integration to compute 
the column densities. 

In our computations we include all reactions listed in the UMIST99 database (Le Teuff 
et al. 2000) for which the reactants and products are species in our set. We adopt the 
recommended UMIST99 rate-coefficients kiji, with several alterations and updates ^. 



^We assume that the abundances of all other elements are negligible in the gas phase, including heavy 
elements such as Mg and Fe which, when abundant, may become important positive-charge carriers in dense 
molecular clouds (Oppenheimer & Dalgarno 1974; dc Boisanger et al. 1996). 

^Our input list of rate coefficients is available at ftp://wise3.tau.ac.il/pub/amicl/pdr. In our compilation 
we use recent determinations of the rate coefficients and branching ratios for dissociative recombinations 
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The depth-dependent FUV photodissociation and photoionization rates are crucial quan- 
tities. They depend on the intensity and spectral shape of the incident radiation field, on 
the photodissociation and photoionization cross-sections, and on the dust grain scattering 
and absorption properties. Dust attenutation is the primary mechanism for the reduction 
of the photorates. Line-shielding is also important for the special cases of H2 and CO pho- 
todissociation ^. For the depth dependent photodissociation and photoionization rates, we 
adopt the biexponential representations 

Fi = dexpi-aiAv - (5iAl) (2) 

and coefficients, Cj, ccj, and calculated by SD95 for the radiative transfer of a Draine field 
penetrating a cloud of large {Aytot = 100) total optical depth (Roberge et al. 1991). We 
assume a dust-to-gas ratio such that Ay = 4.75 x 10~^^(iVH + 2iVH2) (SD95; see also Draine 
2003) where A^h and A^H2 a^'c the atomic and molecular hydrogen column densities (cm^^). 
In Figure 1 we plot the CN and HCN photodissociation rates as functions of Ay, for x = 1 
and a total cloud thickness Aytot = 100. The much more rapid attenuation of the CN 
photodissociation rate compared to HCN is an important feature of the models, and reflects 
the fact that the energetic photons (> 12.4 eV) required to dissociate CN (van Dishoeck 
1987; Thi et al. 2004) are more readily absorbed by the dust grains. 

The cosmic-ray destruction rates in equation (1) include impact ionizations of H, H2 
and He, by the primary cosmic-rays and secondary electrons, and induced FUV photodis- 
sociation and photoionization of the heavier molecules (Gredel et al. 1989). In most of our 
models we adopt a "dense cloud" H2 cosmic-ray ionization rate, ( — 5 x 10~^^ s~^ (e.g., 
Williams et al. 1998), but wc also examine the effect of increasing ( by up to a factor 10 
(McCall et al. 2003; Brittain ct al. 2004). The cosmic-rays are "freely penetrating" and ( is 
independent of cloud depth. Cosmic-ray ionization is the only source of H+, H^, and He+ in 
our models. X-ray ionization (Maloney et al. 1996; Spaans & Meijerink 2005) is excluded. 

The cosmic-ray induced photoionization and photodissociation rates are given by 



(Vejby-Christensen et al. 1997; Larson et al. 1998; Vikor et al. 1999; Jensen et al. 2000; McCall et al. 2003; 
Geppert et al. 2004). For some neutral-neutral reactions we use rate-coefficients listed in UMIST95, or as 
given by Herbst et al. (2000) and Smith ot al. (2004). We do not assume the small (~ 50 K) activation 
barriers postulated by Pineau des Forets et al. (1990) for the neutral-neutral reactions (see §3) involved in 

the nitrogen chemistry. 

•^We set the Doppler parameter equal to 2 km s~^ for all absorption lines. As in SD95, we include H2 
transitions in the Lyman and Werner bands, neglecting the rotational structure. We use the Federman et 
al. (1979) self-shielding formulae for H2, and the van Dishoeck & Black (1989) shielding function for CO. 
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where u is the grain albedo and pi are the photoabsorption "efficiency factors" (Gredel et 
al. 1989). We adopt the efficiency factors Pi hsted in the UMIST99 compilation, and we set 
= 0.5 (SD95). 

Molecule formation requires the presence of H2. We assume that hydrogen molecules 
are formed on grain surfaces with a rate coefficient 

R = 3x 10-^^T^/2^F cm^ s"^ (4) 

where T is the gas temperature, and yp is a "formation efficiency" that depends primarily 
on the grain temperature (Manic et al. 2001; Biham et al. 2001, 2002; Cazaux & Tielens 
2002; Vidali et al. 2004). We set yp — ^ appropriate for low-temperature (~ 10 — 20 K) dust 
grains. We do not consider accretion or ejection of other molecular species onto, or from, 
the grains. 

PDRs consist of several distinct zones or layers with sizes and locations that depend on 
the gas density and FUV field strength. As discussed by SD95, these include the outer H I 
zone, followed by the C II, S II, and Si II zones, terminated by a fully cosmic-ray dominated 
dark core. In the H I zone, rapid photodissociation keeps the hydrogen in atomic form. 
In the C II zone, the hydrogen "self-shields" and becomes molecular, but photoionization 
maintains the carbon as C"*". In the S II, zone the carbon is fully incorporated into CO, but 
photoionization and charge-transfer maintain the sulfur as S"^. In the Si II layer the silicon 
is maintained as Si"*", but the sulfur is either atomic or is incorporated into molecules. The 
relative efficiencies of the various molecular formation sequences, and the resulting molecular 
densities, vary with cloud depth and from zone to zone. For example, in the model considered 
by SD95, for which the hydrogen particle density n = 10^ cm~^ and % = 2 x 10^, the CN/HCN 
density ratio is large (~ 10) in the the H I, C II, and S II layers, but then decreases and 
becomes small (~ 10~^) in the dark core. A large CN/HCN ratio was therefore identified as 
an important feature of (dense) PDRs. SD95 identified additional molecular diagnostic ratios 
of FUV irradiated gas, including OH/H2O, CO+/HCO+, and SO+/SO (see also Sternberg 
et al. 1996). 

The sizes and locations of the various zones also depend on the total gas-phase elemental 
abundances. In our calculations we set the gas-phase elemental abundances equal to the 
interstellar values observed in the diffuse cloud toward ( Oph. Thus, we set C/H=1.32 x 10~^ 
(CardeUi et al. 1993), N/H=7.50 x 10"^ (Meyer et al. 1997), 0/H=2.84 x 10"^ (Meyer et 
al. 1998), Si/H=1.78 x 10"^ (CardeUi et al. 1994), and S/H=8.30 x 10"^ (Lepp et al. 1988). 
We adopt a helium abundance He/H=0.1. In most of our models we keep the gas-phase 
abundances fixed at these values, and independent of cloud depth. In one model we examine 
the effect of reducing the S (and Si) abundances by a factor of 100. 
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In this paper, we study the behavior in cool isothermal clouds. Our results are insensitive 

to the temperature for T < 200 K, and we adopt a representative gas temperature oi T — 
50 K. We consider isobaric clouds, so that the total hydrogen particle densities n = nu + nu.^ 
are constant through the clouds. The densities of all other species increase by a factor of 2 
at the H to H2 transition layer. We present results for n ranging from 10^ to 10^ cm~^. 

Explicit thermal balance computations (Sternberg & Dalgarno 1989; Burton et al. 1990) 
show that in dense (n > lO"' cm~^) PDRs the temperatures in the outer H I zones may 
exceed 10'^ K. Molecular synthesis may then proceed via a "hot gas" chemistry in a thin 
outer layer (e.g. via the efficient production of the reactive intermediate OH). However, the 
contributions to the total CN and HCN column densities from the hot layers are expected 
to be small (SD95). Here we exclude considerations of hot gas, and focus on the chemistry 
in the more extended cooler regions. 

3. CN and HCN Formation and Destruction 

There arc three major pathways to the gas-phase formation of CN and HCN, as illus- 
trated in Figures 2-4. In pathway #1, production occurs via the formation of carbon hydrides 
and intermediate CH and CH2 radicals. In pathway #2 an additional route to CN (but not 
HCN) occurs via oxygen hydrides and the intermediates OH and NO. In pathway #3 the 
molecules are produced via nitrogen hydrides and parent H2CN+ ions. We describe each of 
these pathways in turn. 

3.1. Pathway #1 

In pathway #1, which is important at all cloud depths, CN and HCN are formed via 
the neutral-neutral reactions ^ 

CH + N ^ CN + H (Rl) 

and 

CH2 + N ^ HCN + H . {R2) 



^Reactions (Rl) and (R2) play central roles in the formation of CN and HCN. In our computations we 
assume a rate coefficient of 2.0 x 10^^" cm"^ s^^ at 50 K for (Rl), as implied by pulsed laser experiments 
in room-temperature discharge flows, and a theoretical extrapolation to low temperature (Brownsword et 
al. 1996). This is an order-of-magnitude larger than the rate coefhcient reported by Messing et al. (1981), 
and widely used in many previous PDR models including SD95. For Reaction (R2) we use the theoretical 
value of 5.9 x 10~^^ cm^ s~^ recently calculated by Herbst et al. (2000). 
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The required CH and CH2 intermediates are produced in three ways. First is via radia- 
tive association 

C+ + H2 ^ CH+ + V {R3) 

followed by rapid abstraction 

CH+ + H2 ^ CH+ + H (RA) 

and dissociative recombination 

CH+ + e ^ CH2 + H (i?5) 

CH+ + e ^ CH + H2. (R6) 
Second is cosmic-ray driven proton transfer 

H2 + cr ^ H+ + e (R7) 

H+ + H2 ^ H+ + H {R8) 
C + H+ ^ CH+ + H2 (R9) 

followed by 

CH+ + H2 ^ CH+ + H (RIO) 

and then (R4), (R5), and (R6). 

A third possibility is direct formation ^ via slow radiative association 

C + H2^CH2 + i/ (mi) 

which may become competitive in dense clouds where the fractional ionizations and the 
relative H;^ densities, n^+/n, become small. 

In the outer PDR, the CH and CH2 radicals are removed by FUV photodissociation 
and photoionization. In shielded regions they are removed by 

CH + O ^ CO + H {R12) 

CH2 + O ^ CO + H2 {R13) 
in rapid reactions with the available oxygen atoms. 



^Reaction (Rll) was not included in the SD95 study. Here we adopt fcn = 1.0 x 10 cm^ s ^ for the 
rate coefficient (Smith et al. 2004). 
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3.2. Pathway #2 

Pathway ^^2 is an additional route to CN (but not HCN), and becomes important mainly 
in shielded regions. In this pathway CN is formed via the intermediates OH and NO in the 
reaction pair 

OH + N^NO + H (RU) 

NO + C ^ CN + O . (Rib) 
The main source of OH is cosmic-ray driven proton transfer 

O + H+ ^ 0H+ + H2 (i?16) 

followed by abstractions 

0H+ + H2 ^ H2O+ + H {R17) 

H2O+ + H2 ^ H3O+ + H {R18) 

and dissociative recombination 

H3O+ + e^OH + H + H {R19) 

H3O+ + e ^ OH + H2 . (R20) 

Major removal mechanisms for OH and NO are photodissociation, and reactions with atomic 
oxygen and nitrogen 

OH + O ^ O2 + H {R21) 

NO + N ^ N2 + O {R22) 

in the shielded regions. 

3.3. Pathway #3 

Finally, in pathway #3 both CN and HCN are produced via dissociative recombination 

H2CN+ + e ^ HCN + H {R23) 

H2CN+ + e ^ CN + H2 (i?24) 
where the parent H2CN"'" ions are formed via 

NH + C+ ^ CN+ + H {R25) 

NH2 + C+ ^ HCN+ + H {R26) 
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foUowed by 

CN+ + H2 ^ HCN+ + H {R27) 

HCN+ + H2 ^ H2CN+ + H . (ii:28) 

This pathway is mainly important in shielded regions, but may also operate in the outer H I 
zone. 

In shielded regions, where the N2 density becomes large, NH and NH2 are produced by 
the cosmic-ray driven sequence 

He+ + N2 ^ N+ + N + He (i?29) 

N+ + H2 ^ NH+ + H (i?30) 

NH+ + H2 ^ NH+ + H {R31) 

NH+ + H2 ^ NH+ + H (i?32) 

followed by 

NH+ + e ^ NH + H2 (i?33) 

NH+ + e ^ NH2 + H (i?34) 
In the H I zone, NH may be formed via 

N + H; ^ NH + H (i?35) 

where H2 refers to the sum over all FUV-pumped hydrogen molecules in vibrational states 
with energies large enough to overcome the endothermicity of the reaction. 



3.4. Destruction 

The CN and HCN density profiles also depend on the varying efficiencies of the destruc- 
tion mechanisms. In the outer PDR, photodissociation 

CN + i/ ^ C + N (i?36) 

HCN + ^ CN + H {R37) 
dominates. In shielded layers the CN radical is removed by 

CN -f N ^ N2 + C {R38) 

CN + O ^ CO + N . (i?39) 
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Reaction (R38) dominates at intermediate cloud depths where a high abundance of nitrogen 
atoms are maintained by photodissociation. At large depths (R39) becomes the dominant 
removal process. The saturated HCN molecule is removed primarily by cosmic-ray induced 
photodissociation 

HCN + crp ^ CN + H . (i?40) 

Photodissociation of HCN operates as an important source of CN, both in the PDR and in 
the dark core. Additional removal processes for HCN in the dark core are 

HCN + H+ ^ HCN+ + H (i?41) 

HCN + HCO+ ^ H2CN+ + CO (i?42) 

Reaction (R41) is important in low density clouds (n < 10^ cm~^) where a large proton 
density is maintained by cosmic-ray ionization. Reaction (R42) becomes important in dense 
clouds {n > 10^ cm~^). 



3.5. Free Carbon: C+ and C 

The efficiencies of the CN and HCN formation sequences depend on the availability of 
"free carbon" particles, either as C"*" ions or C atoms. As we discuss in §4, the CN and HCN 
profiles are closely related to the C+ and C density distributions. 

In the outer parts of the PDR, C"*" is produced by FUV photoionization 

C + i/ ^ C+ + e . (i?43) 

In shielded regions the main source is cosmic-ray driven helium-impact dissociative ionization 
of CO 

He + cr ^ He+ + e (RU) 

He+ + CO ^ C+ + O + He . (i?45) 

In the outer PDR, where the electron density is high, the C+ ions are removed by radiative 
recombination 

c+ + e ^ c + 1/ . (me) 

In shielded regions, radiative association (R3), and charge transfer neutralization 

C+ + S ^ S+ + C , (i?47) 

become important. Charge transfer with atomic sulfur plays a particularly important role 
in the partially shielded layers. In the dark core 

C+ + O2 ^ CO + 0+ (i?48) 
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C+ + 02 ^ C0+ + O (i?49) 

dominate. 

Carbon atoms are produced by radiative recombination and charge transfer ([R46] and 
[R47]), FUV photodissociation 

CO + 1/ ^ C + O (i?50) 
and cosmic-ray induced photodissociation 

CO + crp ^ C + O . {R51) 

The C atoms arc removed by photoionization in the outer PDR, by proton transfer and 
radiative association ([R9] and [RH]), and cosmic-ray induced photoionization 



at intermediate depths, and by 



C + crp ^ C+ + e {R52) 

C + O2 ^ CO + O (i?53) 

C + SO ^ S + CO (i?54) 

C + SO ^ CS + O (i?55) 



in the shielded cores. 



4. Reference Model 

The results of our model calculations are displayed in Figures 5, 6, and 7, which show 
cuts through the x, n, ( parameter space we have considered. For each set of cloud parame- 
ters we display the 0+, C, CO, CN, and HCN densities (cm~^) and column densities (cm~^), 
as well as the CN/HCN density and column density ratios, as functions of Ay. The column 
densities are integrated from the cloud surface, in the normal direction. We also display the 
electron density profiles. 

We first discuss our results for a specific "reference model", with n — 10^ cm~^, x — 
2 X 10^, and ( — 5x 10~^^ s~^. We use this model to analyze the depth dependent formation 
sequences and resulting density profiles. We also use it as a point of comparison as the cloud 
parameters are varied. Our reference model is displayed in panels (c), (d), and (e) of Figure 
5, and again in Figures 6 and 7. 

In our reference model, and in our parameter study we adopt T = 50 K. For our 
assumed set of reaction rate coefficients, the results we describe below are quite insensitive 
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to the assumed gas temperature, provided the gas is cold with T < 200 K. For example, 
relative to our results for 50 K, the CN and HCN densities change by less than a factor of 2 
at all if T is decreased to 20 K, or increased to 200 K. 

4.1. C+, C, and CO Density Profiles 

The C+, C, and CO density profiles in our reference model show the familiar structure 
seen in numerous PDR calculations (Tielens & HoUenbach 1985; van Dishoeck & Black 1988; 
Le Bourlot et al. 1993a; Flower et al. 1994; SD95; Jansen et al. 1995; Bakes & Tielens 1998). 
The carbon is ionized in the outer layers, and undergoes a conversion to CO as the ionizing 
FUV radiation is attenuated. An important feature is the double peaked atomic carbon 
density profile. 

The first C peak is associated with the sharp transition of the available gas-phase carbon 
from fully ionized to molecular form near the inner edge of the C II zone. In our reference 
model this occurs at Ay = 2. Due to the rapid production of CO at this location, the 
C density at the peak reaches only ~ 20% of the total gas-phase carbon abundance. The 
transition is, effectively, directly from C^ to CO. After the conversion to CO only trace 
quantities of C and C"*" remain. 

An important property of the C"'"/C/CO transition region is that photodissociation of 
CO continues to dominate the production of the free carbon well beyond the point where 
the conversion to CO is complete. For example, in our reference model, photodissociation 
dominates up to Ay — 2.5, well inside the S II zone. Up to this point, C+ continues to 
be produced by photoionization, and the C''"/C density ratio is set by the balance between 
recombination and photoionization. Thus, in the C II zone, where the production rate of 
atomic carbon via recombination is constant (since Ue ~ "^0+), the C density increases as 
the photoionization rate is attenuated. However, once the conversion to CO is complete the 
atomic carbon density drops as the CO photodissociation rate decreases. This behavior gives 
rise to the first C peak. 

We note that for PAH abundances > 10~^ mutual neutralization, C^ + PAH~ 
C + PAH, can compete with radiative recombination (Lepp et al. 1988) leading to a slight 
shift in the positions of the C+/C/CO transition layers (Bakes & Tielens 1998). We do 
not include the effects of PAHs in our analysis. Similarly, the process of "grain-assisted" 
ion-electron recombination may become important when rieOtr/nag < 1, where (equal to 
1.4 X 10~^^ cm~^ s~^ at 50 K, for C+ ions), and ag < 10^^^ cm^ are the radiative and 
grain- assisted recombination rate coefficients. The efficiency of grain- assisted recombination 
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depends on the grain population size distribution and on the grain charge (Weingartner & 
Draine 2001). In the computations we present here we exclude this process. Including it can 
increase the C densities at the first C peak. We have verified that the CN and HCN density 
profiles we discuss below are not affected significantly. 

At larger depths, e.g. for Ay > 2.5 in our model, the rate of CO photodissociation be- 
comes vanishingly small, and the primary source of free carbon becomes cosmic-ray driven 
helium- impact ionization of CO ([R44] and [R45]). The resulting C^ ions are neutral- 
ized by radiative recombination (R46), and by rapid charge transfer of the C^ ions with 
S atoms (R47). The rapid neutralization keeps the free carbon in atomic form, so that 
nc/nc+ 3> 1. Because the rate of helium impact ionization is independent of cloud depth, 
the C density increases with Ay so long as photoionization dominates the removal of the 
carbon atoms, and as the photoionization rate is attenuated. The second C peak occurs 
at the point where radiative association with H2 (Rll), cosmic-ray induced photoionization 
(R52), and proton transfer reactions with H3 (R9) begin to dominate the removal of the 
carbon atoms. In our reference model this occurs at Ay — 3.7. 

The C density at this location may be estimated analytically, 

2p52C/n+ kgn^+/n + ku 

where kg — 2.0 x 10^ and ku = 1.0 x 10~^^ cm^ s~^ are the rate coefficients for (R9) and 
(Rll), the parameter P52 = 1.0 x 10"^ is the efficiency factor for (R52), and Xhc = 0.1 is the 
He abundance. The first two terms in the denominator becomes small for n > 10^ cm~^, and 
nc is then independent of n. In this limit, ric ~ 0.25 cm~^ for ( = 5 x 10"^'' s~^, in good 
agreement with the numerical value of 0.40 cm~^ at the second C peak. Beyond this point 
the C density decreases, particularly as the O2 and SO densities rise and as the removal 
reactions (R53), (R54) and (R55) become more effective. A sharp drop in the C density 
occurs at Ay ~ 7, and the C density finally stabilizes to a value of 3.6 x 10~^ cm~^ in the 
cosmic-ray dominated dark core. The atomic carbon density in the core depends on n and 
C as discussed further in §5. 

The second C peak marks the point where removal of the C atoms by FUV photoioniza- 
tion becomes inefficient compared to chemical reactions and cosmic-ray induced photoion- 
ization. The minimum between the two atomic carbon peaks occurs near to the point where 
helium impact dissociative ionization of CO replaces photodissociation of CO as the primary 
source of free carbon particles. 



(5) 
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4.2. CN and HCN Profiles 

In our reference model a pronounced peak is present in the CN density profile at Ay = 2, 
at the inner edge of the C II zone, and close to the location of the first C peak. A small 
HCN peak is also present at this location. Throughout the C II zone, the CN and HCN 
molecules are produced primarily via pathway #1, where the formation of the CH and CH2 
intermediates is initiated by radiative association of C"*" with H2 (R3). The CN density peak 
occurs where the C"^ density is still close to its maximum possible value, but where the 
photodissociation rates are diminished. The CN density then decreases with cloud depth as 
the C"*" density drops across the C+/C/CO transition layer. 

Because the carbon is fully ionized in the C II zone, molecule formation is "recombination 
limited" in this part of the cloud. The CN density at the inner edge of the C II zone 
may therefore be expected to scale linearly with the gas density n for a fixed value of the 
"ionization parameter" x/n. This may be seen analytically by writing 

ncR ~ r- feu ~ —^fcnXcn^ (6) 

CH CH 

and 

ncN ~ ~ /chXcXn(- ) n . (7) 

X^ CN ^ CH^ CN 

In these expressions, ki = 2.0 x 10~^° and ^3 = 5.7 x 10~^^ cm"' s^^ are the rate coefficients for 
reactions (Rl) and (R3), and x^'cn — 1-7x10"^ and x^'cn — 5.0x10"^° s~^ are the attenuated 
photodissociation rates at the inner edge of the C II zone. The fraction of recombining of 
CH^ ions that fragment to CH is /ch = 0.43, and Xc = 1.32 x lO^^ and Xn = 7.50 x lO^^ 
are the gas phase abundances of carbon and nitrogen. For n = 10*^ cm~'^ and x = 10^ this 
gives ncN = 5.7 x 10~^ cm~^, in good agreement with our numerically computed value of 
9.0 x 10~^ cm~^ at the peak. The photo-attenuation factors (see equation [2]) depend on 
the visual extinctions at the termination points of the C II zones, and therefore scale with 
x/n. Equation (7) shows that for fixed x/n the CN density at the peak is proportional 
to the product of the carbon and nitrogen abundances and to the cloud density n, and is 
independent of x- These simple scalings are verified by our numerical results in §5. 

The CN density decreases as the supply of free carbon produced by CO photodissociation 
diminishes. However, when the primary source of free carbon switches to helium impact 
ionization, the CN density stops decreasing and a pronounced "plateau" appears in the 
density profile. In our reference model the CN density plateau extends from Ay = 3 to 6. 

The CN density profile fiattens for a combination of reasons. First, while the CH and 
CH2 formation efficiencies decrease as the C+ disappears, this is offset by a more rapid decline 
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in the the CH and CH2 photodestruction rates. Second, the rise in the atomic carbon density 

up to the second C peak leads to an increased efficiency of the proton-transfer sequence 
initiated by (R9) in pathway #1. Third, CN formation via pathway #2 also begins to play 
a significant role as the OH and NO densities increase. 

Furthermore, while HCN (formed via [R2]) is removed by photodissociation, at these 
cloud depths CN is removed by neutral-neutral reactions with N atoms (R38). Photode- 
struction of CN is ineffective in the plateau because of the severe attenuation of the CN 
photodissociation rate (see Figure 1). Importantly, photodissociation of HCN becomes a 
major source of CN in the plateau. The HCN density therefore rises as the FUV field is 
attenuated, but because the atomic nitrogen density remains large throughout, the resulting 
CN density remains insensitive to Ay. 

At these cloud depths molecule formation is "ionization limited" since only a trace 
amount of free-carbon is released from the CO molecules by helium impact ionization. Be- 
cause CN is removed in reactions with a neutral species (N atoms), the CN density is pro- 
portional to the cosmic-ray ionization rate and independent of the gas density. For a simple 
analytic estimate we assume that CN formation is initiated by proton transfer in reactions 
of C with H3 (R9), or by radiative association of C with H2 (Rll), in pathway #1. We also 
assume and that every HCN formation event leads to CN via photodissociation. We assume 
further that the CH and CH2 intermediates are removed by reactions with atomic oxygen 
[(R12) and (R13)], and that these reactions proceed with equal approximate rate coefficients 
ko = 2x 10"^° cm^ s-\ Setting = 1 x 10"^° cm^ s'^ for both (Rl) and (R2) it follows ^ 
that 

ncN ~ 0-5Ci— 1 — TT- (8) 

^38 Ko Xq 

where = 3 x 10"^*^ cm^ s~^ is the rate coefficient for removal reactions of CN with N. 
Evaluating for X^e = 0.1, Xq = 2.84 x 10~^ and C = 5 x 10"^^ s"^ yields ncN = 9 x 10"^ 
cm^"^, in good agreement with the numerically computed ticn ~ 2 x 10^^ cm^"^ in the CN 
plateau (see Figure 5). Thus, ncN is proportional to ^ and independent of n, and with the 
above assumptions ^, is independent of both the gas-phase nitrogen and carbon abundances. 

The CN density plateau is maintained up to the point where photodissociation by the 
incident FUV photons is the dominant HCN removal mechanism, or up to Ay = 6 in our 



^The factor of 0.5 in this expression enters because the cosmic ray ionization rate of He is half that of H2. 

^Thcsc assumptions break down if, for example, Xc > Xq so that C atoms become abundant rather 
than being mainly locked in CO, or if Xq becomes so small such that formation via pathway #1 becomes 
negligible. 
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reference model. The HCN density reaches a peak value of 8.6 x 10~^ cm~^ at this point. At 
larger depths cosmic-ray induced photodissociation becomes the dominant HCN destruction 
mechanism, proceeding at a rate independent of Ay- CN is removed by reactions with 
oxygen atoms (R39). The CN and HCN formation efficiencies both decline as less C and C"*" 
is available, and the CN and HCN densities decrease to their dark core values of 1.2 x 10~^ 
and 2.6 x 10~^ cm~^ respectively. 

Figure 5e shows that the CN column density rise sharply to A^cn = 6 x 10^^ cm^^ at 
the CN peak at Ay = 2, and that much of the CN column is built up at this location. The 
HCN column remains small at Ay = 2, but continues to increase with cloud depth. A^hcn 
approaches 7 x 10^^ cm~^ at Ay = 6 where the HCN volume density is at its maximum. The 
ratio A^cn/-Ahcn ~ 10 at Ay = 2, and decreases to ~ 0.1 at Ay — 6 (see Figure 5f). 

5. Pcirameter Study 
5.1. Density 

Figure 5 displays our results for n equal to 10^, 10^, and 10^ cm~^, in addition to our 
10^ cm~^ reference model, and illustrates the effects of varying the cloud density. In this 
sequence we keep the ionization parameter x/n constant at 0.2 cm~^, so that the C+/C/CO 
transition layers occur at the same location {Ay — 2) in all four models. 

For fixed x/n, the depth at which helium impact ionization of CO replaces photodisso- 
ciation of CO as the source of free carbon, increases with the cloud density n. However, the 
depth at which removal of atomic carbon by photoionization becomes ineffective compared 
to proton transfer or radiative association ([R9] or [Rll]), is independent of n. Therefore, 
the relative height of the second C peak (at Ay = 4.2) decreases with n, and the peak dis- 
appears at sufficiently high densities (> 10^ cm^^). For densities above 10^ cm"'^, riQ ^ 0.4 
cm~^ at this location, consistent with Equation (5). At lower densities all of the available 
carbon is atomic at the location of the second C peak, and ric is then limited by the carbon 
abundance Xq (see Figure 5a). 

At low n, the free carbon density remains relatively high at large Ay and in the dark 
core. At high n the free carbon density becomes small. As we discuss below, this behavior 
reflects the transition from the "high-ionization phase" to "low-ionization phase" in the 
cosmic-ray dominated core (Le Bourlot et al. 1993b; 1995; Lee et al. 2000). The asymptotic 
values of the C densities in the dark core influence the shapes of the atomic carbon density 
proflles at intermediate depths in the PDRs. 
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Figure 5 shows that the CN density at Ay — 2 increases hnearly with n, from ~ 10~^ to 
10"^ cm~^, in accordance with Equation (6). The "CN peak" becomes more pronounced as 
n increases. In the n = 10^ cm~^ model, the density of FUV-pumped H2 becomes large, and 
CN is also formed via pathway ^3 initiated by (R35). For n > 10^ cm~^ the "CN plateau" 
from Ay = 3 to ~ 6 is apparent. At these depths ncN ~ 2 x 10~^ cm~^ independent of n, 
consistent with Equation (8). At low n, the CN density remains large to high Ay due to the 
elevated free carbon densities. 

For fixed x/""-) FUV photodissociation of HCN is effective to greater depths as n in- 
creases. The location of the HCN density peak therefore moves from a visual extinction of 
5.6 to 7.2 as the gas density is increased from 10^ to 10^ cm~^. The column density ratio 
A^cn/A^hcn remains large to greater depths as n is increased, but the density ratio ncN/^HCN 
falls off sharply with Ay (see Figure 5). 

5.2. FUV Intensity, Sulfur Abundance, and Cosmic-Ray Ionization Rate 

In Figure 6, we set n equal to 10^ cm~^, and we vary x from 2 x 10^ to 2 x 10®. 
This illustrates the effects of varying the incident FUV field intensity. The location of the 
C+/C/CO transition layer moves from Ay — 1.3 to 3.5 for this density and range of FUV 
intensities. The position of the second C peak is also shifted to larger Ay, as is the location 
of the sharp drop in the C density, which moves from Ay = 6 to 9. The approach to the 
dark-core conditions occurs at greater Ay as x is increased. 

Correspondingly, the locations of the CN peak and the inner HCN peak also move 
inward as x is increased. Because the CN photodissociation rate declines most rapidly with 
Ay (see Figure 1) ncN/n-ncN and A^cn/-/Vhcn at the CN peak increase with x, from ~ 10 to 
40 in our models. 

In Figure 7, we show the effects of a reduced gas-phase S (and Si) abundance, and an 
increased cosmic-ray ionization rate, for our n — 10^ cm~^ and x — 10^ model. In panels (a), 
(b), and (c) of Figure 7, the S (and Si) abundances are reduced by a factor of 100. Charge 
transfer (R47) between C"*" and S is then less effective as a neutralizing mechanism, and the 
second C peak at large Ay is reduced. The formation efficiency of HCN in pathway 7^1, via 
radiative association of C with H2 (R-H) followed by the reaction of CH2 with N (R2), is 
therefore also diminished, leading to a reduction of the inner HCN density peak (compare 
Fig. 7a and Fig. 7d). The resulting CN/HCN ratios are therefore enhanced in the CN plateau 
region, and in the dark core. 

The bottom two rows of Figure 7 show the effects of increasing the cosmic-ray ionization 
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rate, by factors of 3 and 10, to 1.5 x 10~^^ and 5 x 10~^^ s~^. Increasing ( mainly affects 
the atomic carbon density profile. The first C peak is unaltered because it is is controlled 
by purely photo-processes. However, the second C peak is enhanced by the increased rate 
of helium impact destruction of CO (R44). The behavior is consistent with Equation (5). 
Initially, the C density increases linearly with (. However, for sufficiently large ( the atomic 
carbon is removed by cosmic-ray induced photoionization (R52) and proton transfer reactions 
with Hg (R9), rather than by radiative association with H2 (Rll). In this limit the C density 
is independent of (, and all of the available gas phase carbon becomes atomic, as occurs for 
our ( — 5 X 10"^^ model. 

The free carbon densities remain large throughout the cloud for high (. As we now 
discuss, the shift from a low to high free carbon density with increasing ( is related to a 
"phase change" that occurs in the dark core. 

5.3. Dark Core 

In Figure 8 we display the densities, rii, and density fractions, rii/n, for C+, C, CO, CN, 
HCN, and free electrons, for fully opaque dark core conditions. Externally incident FUV 
photons are excluded, and the chemistry is driven entirely by cosmic-ray ionization. For 
these conditions, it follows from Equation (1) that the density fractions rii/n depend on a 
single parameter, the ratio of the cloud density to the ionization rate. In Figure 8 we plot 
the solutions as functions of n/C_i7, where C-17 is the cosmic-ray ionization rate normalized 
to 1.0 X lO"^'^ s-^ The phenomenon of "bistability" (Le Bourlot et al. 1993b, 1995; Lee 
et al. 2000) is apparent in Figure 8, where for densities between 370 and 675 cm^'^ two 
stable solutions exist. A high-ionization phase (HIP) occurs at low n, and a low-ionization 
phase (LIP) occurs at high n. The two phases may coexist where the gas is bistable. We 
will present our own discussion of bistabihty elsewhere (Boger & Sternberg 2005). Because 
bistability occurs for a narrow range of densities the transition from the HIP to LIP may 
be said to occur near a "critical density" ncrit/C-i7 ~ 500 cm~^ (for our assumed gas-phase 
abundances). In our computations, Ue/n > 10~^ in the HIP, and Ue/n < 10~^ in the LIP, 
and the fractional ionization drops by a factor ~ 10, at the transition density. The C"*" and C 
densities are large in the HIP, with ric/nco ~ 1) and small in the LIP, with ric/nco < 10~^, 
consistent with the findings of Flower et al. (1994). Correspondingly, riQ-^/n and nncN/''^ are 
large for n < ricnt, and small for n > Hdit- Furthermore, ncN/'^HCN is large in the low-density 
HIP, and becomes small in the high-density LIP. In the HIP, ncN/^HCN decreases from ~ 20 
to 0.07 for n/C-17 ranging from 10 to 675 cm""^. In the LIP, ncN/^HCN decreases from 0.4 
to 3 X 10~^ cm~^ for n/C-17 ranging from 370 to 10^ cm~^. 
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Returning now to Figure 7, we note that the increase in the free carbon density at large 

Ay that occurs as ( is increased, is due to the rise in ncnt, and a corresponding transition 
from LIP to HIP conditions, for the specific gas density of 10*^ cm^^. The phase transition 
is affected by the presence of some FUV radiation. For example, for ( = 1.5 x 10~^^ s~^ it 
occurs at Ay ~ 7 rather than in the fully opaque core (see Figure 7g). As the ionization 
rate is increased further, the gas is converted to the HIP throughout (see Figure 7j). 

6. Discussion and Summary 

CN and HCN molecules have been observed in a diversity of Galactic and extragalac- 
tic sources, and the CN/HCN intensity and density ratios have been used as diagnostic 
probes of FUV irradiated molecular gas. One recent and interesting example is the well- 
known starburst galaxy M82 {d = 3.9 Mpc; L = 3.7 x 10^° L©; e.g. Forster-Schreiber et 
al. 2003). Fuente et al. (2005) have reported measurements of CN 1-0 (113.490 GHz), 
CN 2-1 (226.874 GHz), and HCN 1-0 (88.631 GHz) hne emissions across the inner 650 pc 
star-forming molecular disk in M82. They find that A^cn/-/Vhcn ~ 5, and argue that the large 
ratio is indicative of a giant and dense PDR bathed in the intense field of the starburst. They 
derive n ~ 10^ to 10^ cm~^, and x ~ 10^- Fuente et al. (2005) conclude that Ay < 5 in 
the M82 clouds, since for optically thicker clouds the CN/HCN column density ratio would 
be smaller than observed. Our results provide support for these conclusions. In this picture 
~ 10 to 20 individual clouds along the line of sight are required, since for characteristic 
elemental abundances the computed CN and HCN columns for Ay < 5 (see Figures 5 and 
6) are about an order of magnitude smaller than observed (A^cn = [2 ± 0.5] x 10^^ cm~^ and 
A^HCN = [4 ± 0.5] X 10^3 cm-2). 

Alternatively, the large CN/HCN ratio might be a signature of a very large cosmic-ray 
ionization rate, up to ~ 5 X 10-^^ s-^ as invoked by Suchkov et al. (1993) for M82 and 
other starburst galaxies. For such high ionization rates the CN/HCN ratio could remain > 1 
even in dense and opaque cores if these are maintained in the high-ionization phase. Our 
results indicate that for ^ = 5 x 10~^^ s~^, the HIP is maintained up to n ~ 4 x 10^ cm~^ 
(see Figure 8). This possibility is perhaps more compatible with the large C/CO ratio 
> 0.5 inferred from observations of the 492 and 809 GHz C I fine-structure lines in M82 
(Schilke et al. 1993; Stutzki et al. 1997; see also Gerin & Phillips 2000). As discussed above, 
a large C/CO ratio is a signature of the HIP. 

In this paper we have presented a theoretical study of CN and HCN molecule formation 
in dense interstellar clouds exposed to intense FUV radiation fields. We have analyzed the 
behavior of the CN/HCN density ratio for a wide range of conditions, with the aim of showing 
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how this molecular ratio may be used as a diagnostic probe of molecule formation in FUV 
irradiated gas. For this purpose, we have constructed detailed models in which we solve the 
equations of chemical equilibrium as functions of optical depth, for uniform density clouds 
at constant gas temperature (50 K). Our results are insensitive to the gas temperature for 
cold clouds with T < 200 K. We consider clouds with hydrogen particle densities ranging 
from n = 10^ to 10^ cm~^, and FUV radiation intensities ranging from % = 20 to 2 x 10^, 
appropriate for star-forming clouds near young OB stars and clusters. We present results for 
cosmic-ray ionization rates ranging from 5 x 10~^^ to 5 x 10~^^ s~^, and we also examine the 
effects of large (factor 100) sulfur depletions on the computed density profiles. We present 
calculations of the density profiles for CN and HON, and for the associated species C+, C, 
and CO. We analyze the behavior from the outer FUV photon-dominated regions, into the 
fully opaque cosmic-ray dominated cores. 

In this paper, we adopt a fixed characteristic shape for the FUV radiation spectrum, 
and we do not examine the possible signatures of "soft" versus "hard" FUV fields (van Zadel- 
hoff ct al. 2003). Time-dependent effects due to episodic shadowing (Storzer et al. 1997), 
fluctuating radiation fields (Parravano et al. 2003), grain processing (Charnley et al. 2001), 
turbulent mixing and diffusion (Papadopoulos et al. 2004), and other dynamical processes, 
are not considered here. 

Our models show how observations of the CN/HCN abundances ratio in molecular 
clouds may used as probes of FUV and cosmic-ray driven gas-phase chemistry for a wide 
range of conditions. We find that in dense gas, CN molecules arc characteristically and 
preferentially produced near the inner edges of the C II zones in the PDRs. This is where C"*" 
begins to recombine and where atomic carbon is incorporated into CO. Molecule formation is 
"recombination limited" at these depths, and for fixed x/^? the CN density is proportional to 
the cloud density, and to the gas-phase carbon and nitrogen abundances. For n>10^ cm~^, 
and for clouds with linear sizes corresponding to visual extinctions Ay < 10, the entire 
integrated CN column density is built up at the C^/C/CO transition layer. For characteristic 
interstellar carbon and nitrogen gas-phase abundances the predicted CN columns are ~ 
3 X 10^^ cm~^. HCN is rapidly photodissociated in the outer parts of the PDRs including 
the C"'"/C/CO transition layers. Because HCN is more vulnerable to photodissociation, the 
CN/HCN density ratio is large at low Ay, and decreases with increasing optical depth. 
We find that the CN/HCN density ratio typically decreases from > 10 in the C+/C/CO 
transition layers to < 0.1 in the opaque cores. 

At intermediate and large depths, the bulk of the gas-phase carbon is locked in CO 
molecules, and the CN and HCN densities depend on the rate at which carbon is released from 
CO by cosmic-ray driven helium impact ionization. Molecule formation is then "ionization 
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limited" and occurs with an efficiency proportional to the cosmic-ray ionization rate. In dense 
clouds, an enhanced abundance of atomic carbon is maintained at intermediate depths, where 
charge transfer of C"*" with S is effective, and where photo-destruction reduces the efficiency 
with which the carbon atoms are removed in reactions with molecules. At these intermediate 
depths the CN density is insensitive to Ay- However, the HCN densities increase with Ay 
as the destructive FUV photons are absorbed. 

The C^, C, CO, CN, and HCN densities in the opaque cores depend on whether the 
cores are in the "low-ionization" or "high-ionization" phases that arc possible for such gas. 
We present computations for gas in both phases. The transition from the HIP to LIP occurs 
at a critical value n/Q-n = 10^ cm~^, consistent with previous findings. The CN/HCN 
density ratio can become large ^ 1 in the LIP, but remains small in the HIP. 

For dense gas in Galactic molecular clouds, and for characteristic cosmic-ray ionization 
rates, a large CN/HCN density ratio may be interpreted reliably as an indicator of molecule 
formation in PDRs. In clouds exposed to enhanced fluxes of cosmic-rays, as perhaps occurs in 
starburst galaxies, a high CN/HCN ratio may alternatively be an indicator of opaque clouds 
in the high-ionization phase. Measurements of the C/CO ratio can be used to distinguish 
between these two possibilities. 
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Fig. 1. — Photodissociation rates of CN and HCN as functions of visual extinction for a unit 
Draine field. 
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Fig. 2. — CN and HCN formation pathway #1 via carbon hydride intermediates. 
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Fig. 5. — Model results for n equal to 10^, 10^, 10^, and 10^ cm~^ (from top to bottom) 
for constant ionization parameter x/^ = 0-2 cm^. The cosmic-ray ionization rate ( = 
5 X 10~^^ s~^ in all four models. Displayed profiles are for C^, C, CO, CN, HCN, and 
electrons, as functions of Ay. The lefthand panels display the volume densities. The middle 
panels show the integrated column densities. The righthand panels show the volume (dashed) 
and column (solid) CN/HCN density ratios. The second row, with n = 10^ cm~^, % = 2 x 10^ 
is our "reference model" . 
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Fig. 7. — Results for n = 10^ cm ^, and x = 2 x 10^. In the upper row the gas phase 
abundances of S and Si are reduced by a factor of 100. The second row is our reference 



model, with standard abundances and C = 5 x 10 s . In the third and fourth rows, ( is 
increased to 1.5 x 10~^^ and 5 x 10~^^ s"^. 
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Fig. 8. — Densities of CN, HCN, C^, C, CO, and electrons, in the cosmic-ray dominated 
cores, showing the high and low ionization phases. Panels (a) and (b) show the densities 
ni/C-17 and density fractions rii/n, as functions of n/C-17, where C-17 is the cosmic-ray 



